Noise Thermometry with Two Weakly Coupled Bose-Einstein Condensates 
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Here we report on the experimental investigation of thermally induced fluctuations of the relative 
phase between two Bose-Einstein condensates which are coupled via tunneling. The experimental 
control over the coupling strength and the temperature of the thermal background allows for the 
quantitative analysis of the phase fluctuations. Furthermore, we demonstrate the application of these 
measurements for thermometry in a regime where standard methods fail. With this we confirm that 
the heat capacity of an ideal Bose gas deviates from that of a classical gas as predicted by the third 
law of thermodynamics. 

PACS numbers: 05.30.Jp, 03.75.-b, 05.40.-a, 74.40.+k 
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The generation of two independent matter-wave pack- 
ets by splitting a single Bose-Einstein condensate (BEC) 
is a well established technique 0, 0, El in the field of 
atom optics. New phenomena arise if the two separated 
parts can still coherently interact in analogy to Joseph- 
son junctions in condensed matter physics and su- 
perfluid Helium Josephson weak links |sjj. An advantage 
of the realization of weakly coupled BEC in a double- 
well potential is the possibility to observe the phase 
difference between the two macroscopic wave functions 
directly. Our experimental investigation of this relative 
phase reveals that it is not locked to zero but exhibits 
fluctuations. Two fundamental types of fluctuations are 
discussed in the literature, quantum fluctuations Q and 
thermally induced fluctuations §. In this letter we re- 
port on the experimental investigation of thermal fluctu- 
ations of the relative phase arising from the interaction of 
the BEC with its thermal environment, which is always 
present. 

The essential prerequisite for the investigation of these 
thermally induced phase fluctuations is the ability to pre- 
pare a BEC adiabatically in a symmetric double- well po- 
tential and to adjust its temperature. In our experiments 
this is achieved by splitting a single 87 Rb BEC produced 
and trapped in an optical dipole trap by slowly ramp- 
ing up a barrier in the center. The tunneling coupling 
is adjusted by the barrier height and its strength can be 
deduced from numerical simulations of the BEC in the 
trap using the model described in 0. The temperature 
of the BEC is adjusted by holding the cloud in the trap, 
where due to fluctuations of the trap parameters energy 
is transferred to the atoms. Once the final temperature 
is reached a standing light wave is ramped up generating 
a barrier in the center, leading to an effective double-well 
trapping potential (upper part of Fig. la). 

When the potential is switched off the matter-wave 
packets start to expand, overlap and form a double- 
slit interference pattern which depends on their relative 
phase as indicated in the lower part of Fig. la. Repeat- 
ing the interference measurements reveals that this rel- 
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FIG. 1: Observation of thermal phase fluctuations. The ex- 
perimental steps are depicted in (a). A Bose-Einstein conden- 
sate (solid line) is prepared in a double- well potential (dashed 
line) by adiabatically ramping up the barrier. The relative 
phase can be measured after a time-of-flight expansion by an- 
alyzing the resulting double-slit interference patterns (black 
line) . (b) shows polar plots of the relative phase obtained by 
repeating the experiment up to 60 times. The graphs show 
measurements for four different temperatures T at constant 
tunneling coupling energy Ej, i.e. constant barrier height. 
The phase fluctuations increase with increasing temperature, 
(c) shows polar plots of the relative phase for a constant tem- 
perature at four different tunneling coupling energies, i.e. dif- 
ferent barrier heights. Here the fluctuations are reduced with 
increasing coherent tunneling coupling showing the stabiliza- 
tion. 



ative phase is not constant but fluctuates around zero. 
The general behavior of these phase fluctuations is con- 
nected to two parameters: the temperature of the sys- 
tem randomizing the phase and the tunneling coupling of 
the two matter-wave packets stabilizing the phase. The 
results depicted in Fig. lb show that the phase fluctu- 
ations become more pronounced as the temperature is 
increased since the fluctuations outweigh the stabilizing 
effects. From this point of view it is expected - and also 
experimentally observed (Fig. lc) - that keeping the tern- 
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perature constant and increasing the tunneling coupling 
leads to a reduction of the fluctuations. A measure for 
the fluctuations is the coherence factor a = (cos^) Q 
which is directly connected to the visibility of the ensem- 
ble averaged interference fringes. 
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FIG. 2: Loss of the coherence of the bosonic Josephson junc- 
tion due to the coupling to a thermal environment, (a) shows 
the transition from coherent single realizations to incoherent 
ensemble averages. In the single realization a clear interfer- 
ence signal is observed, where the visibility is decreased due to 
the finite optical resolution of the imaging setup. After aver- 
aging over 10 realizations the visibility is reduced and after av- 
eraging over 50 realizations the coherence is lost (a = 0.046). 
In (b) the coherent evolution of the bosonic Josephson junc- 
tion is depicted. The results shown are obtained by repeating 
the experiment at the same temperature as above but at a 
stronger tunneling coupling. Here the averaging leads to only 
a small degradation of the visibility (a = 0.87) showing how 
coherent coupling can counteract dephasing processes. 

The coherence of the system can be visualized as shown 
in Fig. 2. For every single realization below the critical 
temperature the experiment reveals interference patterns 
with high visibility. However, the visibility is reduced by 
averaging over many realization and for high tempera- 
tures it disappears completely as the mean fluctuations 
of the relative phase become comparable to it. The loss of 
coherence due to thermal fluctuations is shown in Fig. 2a 
corresponding to a = 0.046. At the same temperature 
the coherence of the weakly coupled condensates can be 
maintained by increasing the coupling. For a tunneling 
coupling energy larger than the thermal energy the phase 
is locked to zero and the averaging reduces the visibility 
only slightly as shown in Fig. 2b where the coherence 
factor is given by a = 0.87. The dependence of the co- 
herence factor on the two parameters is the consequence 
of a universal scaling law which can be explained by the 
classical model discussed in the following. 

The dynamics can be described in terms of a two 
mode approximation by assuming weak coupling between 
the localized modes of the BEC. This corresponds to 
treating the BEC in the double-well potential as two 
separated matter-wave packets connected via tunneling 



through the barrier. In the following we will use the 
acronym for bosonic Josephson junction (BJJ) to de- 
scribe this system. Within the two mode approximation 
the dynamics of the BJJ can be described by two conju- 
gate variables, the atom number difference between the 
matter-wave packet on the left (1) and on the right (r) 
An = (Ni — N r )/2 and their relative phase cf> = <p r — <fii 
111 LLil lilt Il2j | . The Hamiltonian governing the evolution 
of the two conjugate variables in the limit of small An is 
given by 



H = ^An 2 
2 



E* ■ cos < 



(1) 



where E c accounts for the atom-atom interaction in both 
condensates and Ej is the tunneling coupling energy re- 
sulting from the spatial overlap of the wave functions. 
This Hamiltonian also describes the classical motion of a 
particle with mass 1 / E c and momentum An at position <f> 
in a periodic potential. In our experiments with temper- 
atures T > lOnK the quantum fluctuations 0, [13,13113 
are small compared to the thermal fluctuations and there- 
fore are neglected. Their influence can be estimated in 
the limit of small <p m which Eq. (1) can be approximated 
by a harmonic oscillator with the characteristic quantum 
mechanical energy splitting hu> p = y/E c ■ Ej where uj p is 
the plasma frequency, leading to the quantum mechan- 
ical fluctuations of both variables: (An 2 ) w */ Ej/4E C 
and (<P 2 ) « y/Ec/AEj. 

The system variables Ej and E c can be calculated 
from the experimental parameters. The trapping fre- 
quencies of the three-dimensional harmonic trap are lo x = 
2tt ■ 90(2)Hz and w VtZ = 2tt • 100(2)Hz. The periodic 
potential of V = Vb/2(1 + cos(27r/A • x)) is realized by 
the interference of two laser beams at a wavelength of 
830nm crossing under an angle of 10° resulting in a stand- 
ing light wave with periodicity of A = 4.8(2)/im and 
is ramped up to a height of Vo/h = 500Hz to 2500Hz. 
The number of atoms in the BEC fraction is chosen to 
be 2500(500). After the preparation of the BEC in the 
double-well trap the relative phase of the two matter- 
wave packets is measured by analyzing the double-slit 
interference patterns formed after time-of-flight of 5 and 
6ms. The visibility of these patterns is reduced due to the 
short expansion time and the finite optical resolution of 
the imaging system. Further details of the experimental 
setup can be found in . 

The relevant quantities can be calculated from these 
parameters using the improved two mode model 
E c /kB is on the order of 20pK and Ej/ks is between 
30pK and 400nK leading to ftw p jkB being between 
25pK and 3nK. Thus, both necessary conditions for the 
classical limit are fulfilled: Ej E c leading to small 
quantum fluctuations of <j> and E c 3> Ej/N 2 (where N is 
the total number of atoms in the BEC) leading to small 
quantum fluctuations of An/N. Hence, our experiment 
can be discussed in the classical framework where the 
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thermally induced phase fluctuations are closely analo- 
gous to the Brownian motion of a particle in a sinusoidal 
potential. 




FIG. 3: Scaling behavior of the coherence factor of a bosonic 
Josephson junction. Each point is obtained by averaging over 
the cosine of the phases of at least 28 (in average 40) mea- 
surements at the same experimental conditions. The coher- 
ence factor a is plotted as a function of the scaling param- 
eter ksT/Ej which is varied over three orders of magnitude 
(49nK < T < 80nK, 0.6nK < Ej/k B < 300nK). It shows 
good agreement with the theoretical prediction of the classi- 
cal model Eq. (2) indicated by the solid line where also the 
uncertainty arising from the fitting error of the phase is taken 
into account. Typical error bars are shown which result from 
statistical errors and uncertainties of the experimental param- 
eters (potential parameters, atom numbers, temperature). 

For a quantitative analysis in the thermodynamic limit 
at UbT 3> huip the coherence factor |8( can be calculated 
by a thermal average assuming a Boltzmann distribution 
for the relative phases 

d(f) ■ cos (f> ■ exp(Ej/ksT ■ cos </>) 

a=<coscf)>= — - . 

J d<j) ■ exp(Ej/kBT ■ cos 4>) 

(2) 

Eq. (2) points out that the relevant scaling parame- 
ter for thermal fluctuations is the ratio between thermal 
energy fc^T and tunneling coupling energy Ej. Fig. 3 
shows the experimentally obtained coherence factors as 
a function of this scaling parameter. Every data point 
represents on average 40 measurements. In these experi- 
ments the temperature of the system is changed between 
49nK and 80nK by evaporatively cooling the sample to 
the lowest temperature and subsequently increasing the 
temperature by holding the atoms in the trap for differ- 
ent times. The temperature of the sample is measured 
with the standard time-of-flight expansion method. The 
tunneling coupling energy is varied between 0.6nK-fcs 
and 300nK-fcs by adjusting the height of the potential 
barrier. Ej is obtained from numerical calculations us- 
ing independently measured trap parameters and atom 
numbers. It is important to note that the recently de- 
veloped improved two mode model is used for these 



calculations because it leads to quantitative agreement 
between theoretical predictions and experimental mea- 
surements of dynamical quantities The solid line 
corresponds to the theoretical prediction of the classi- 
cal model (Eq. (2)) where all parameters are determined 
independently. It also includes the fitting error of the 
relative phase which arises from the finite optical resolu- 
tion and leads to a reduction of the coherence factor. As 
shown in Fig. 3 the general behavior of the coherence is 
confirmed over a three orders of magnitude variation of 
ksT/Ej. These measurements reveal that the BJJ has a 
higher degree of coherence than expected. This deviation 
might possibly be explained by an increase of the tunnel- 
ing coupling resulting from the excitation of transverse 
modes with higher energies which are neglected by the 
two mode approximation. 

Independent measurements have been performed for 
the lowest temperatures (T = 15nK) to test for thermal 
equilibration. The measurements of a were compared for 
different Ej for two ramping schemes. The first scheme 
was ramping up the barrier in 1.3s and the second scheme 
was holding the atoms for Is in the trap and then ramp- 
ing up the barrier within 0.3s. For Ej/kB > InK both 
schemes lead within the experimental errors to the same 
results. Thus for the fluctuation measurements the ramp- 
ing in 300ms is expected to be adiabatic with respect to 
the response time of the BJJ given by the inverse plasma 
frequency and thus ensures the thermal equilibrium. 

In the following we present the application of the phase 
fluctuation measurements for thermometry far below the 
critical temperature of Bose-Einstein condensation (T c ). 
The temperature of the system can be directly deduced 
from the variance of the phase if the tunneling coupling 
is known. In order to apply the phase fluctuation mea- 
surements for thermometry we introduce an empirical 
effective tunneling coupling Ej to account for effects 
beyond the classical approach. For the range of 25nK 
< Ej/kB < 90nK we deduce from the results shown 
in Fig. 3 that Ef f = 1.33 • Ej. The fundamental dif- 
ference between this method and previous suggestions 
using phase fluctuations of elongated Bose-Einstein con- 
densates for thermometry is that the BJJ is not re- 
stricted to a quasi one dimensional situation but can be 
employed for all geometries. Furthermore, this method 
can be applied for all temperature ranges by tuning E c 
and Ej such that thermal effects dominate and quantum 
fluctuations are negligible. 

As a proof of applicability of this new type of ther- 
mometer we observe how the temperature of a BEC in 
a harmonic trap increases in time (see Fig. 4), which re- 
veals clearly the effect of quantum statistics below the 
critical temperature. In these experiments the lowest 
temperatures (T < T c /3) can only be measured with the 
phase fluctuation method since the thermal fraction is 
too small to be observed in time-of-flight measurements 
(less than 100 atoms with about 2500 atoms in the BEC 
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FIG. 4: Heating up of a Bose gas. The filled circles cor- 
respond to measurements employing the phase fluctuation 
method and the open circles to the results obtained with 
the standard time-of-flight method. The grey shaded region 
shows the critical temperature expected from the experimen- 
tal parameters and their uncertainties. The solid line is a 
fitting function assuming a power law for the heat capacity 
C cc (T/T c ) d of the Bose gas below the critical temperature 
T c — 59(4)nK, a constant heat capacity above the critical 
temperature and a temperature independent transfer rate of 
energy. From this fit we deduce d = 2.7(6) which is con- 
sistent with the theoretical prediction of d — 3 for an ideal 
Bose gas in a three-dimensional harmonic trap. The dashed 
line represents the expected behavior of an ideal classical gas 
for increasing temperature which makes the difference arising 
from quantum statistics evident. 



fraction). For longer heating times the standard time- 
of-flight method can be applied and confirms the consis- 
tency of the two approaches in the overlap region. The 
solid line corresponds to a fitting function for the tem- 
perature where we assume a mean critical temperature of 
T c = 59nK (deduced from independent measurements), 
a temperature independent transfer rate of energy per 
particle and a power law for the temperature dependent 
heat capacity C = (d + 1) • C t h ■ {T/T c ) d where C th is 
the heat capacity of a classical gas. The shown excellent 
agreement is obtained for a heating rate of 2.3(2)nK/s 
for a classical gas and d = 2.7(6). Thus the expected ex- 
ponent d = 3 for an ideal Bose gas in a three-dimensional 
harmonic trap [lflj is experimentally confirmed. The ex- 
pected increase of temperature of a classical gas is indi- 
cated by the dotted line and shows clearly the difference 
between the quantum and the classical behavior of ideal 
gases. 

In summary, we have presented a quantitative analy- 
sis of thermally induced phase fluctuations in a bosonic 
Josephson junction. Our observations show that a uni- 
versal scaling law describes the behavior of the coherence 



and its control leads to new applications. A method is 
presented for ultra-low temperature measurements, with 
which we have confirmed that the heat capacity of a de- 
generate Bose gas vanishes in the zero temperature limit 
as predicted by the third law of thermodynamics [2(J . 
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